Unextendible product bases and locally unconvertible bound entangled states 
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Mutual convertibility of bound entangled states under local quantum operations and classical 
communication (LOCC) is studied. We focus on states associated with unextendible product bases 
(UPB) in a system of three qubits. A complete classification of such UPBs is suggested. We prove 
that for any pair of UPBs S and T the associated bound entangled states ps and pr can not be 
converted to each other by LOCC, unless S and T coincide up to local unitaries. More specifically, 
there exists a finite precision e(S, T) > such that for any LOCC protocol mapping ps into a 
probabilistic ensemble (p a ,p a ), the fidelity between pr and any possible final state p a satisfies 
F(p T , Pa )<l-e(S,T). 
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INTRODUCTION AND SUMMARY OF 
RESULTS 



One of the most challenging problems in the field of 
quantum information is to understand what transfor- 
mations of multipartite entangled states can be realized 
by local quantum operations and classical communica- 
tion (LOCC). This problem has many important appli- 
cations such as construction of entanglement distillation 
protocols quantum data hiding [2(, and separability 
check 3]. The following question is of particular im- 
portance for all these applications. Given mixed states 
ps and pT of some fixed composite quantum system, is 
it possible to convert ps into pr by LOCC transforma- 
tions? More specifically, we will be interested in stochas- 
tic approximate convertibility. Recall, that an output 
of a LOCC transformation applied to the state ps is a 
probabilistic ensemble (p a , p a ), where a is a classical ran- 
dom variable describing results of intermediate measure- 
ments and having a probability distribution {p a }- The 
state p a is an output quantum state conditioned on a. 
By stochastic convertibility we mean that p a must coin- 
cide with pt for at least one outcome a (we assume that 
p a > for all a). In other words, the conversion suc- 
ceeds with some probability which is generally smaller 
than one. By stochastic approximate convertibility we 
mean that for any given precision e > the state ps can 
be stochastically converted into some state p £ satisfying 
F(pT,p e ) > 1 — e, where F is the fidelity. Note that 
the probability for the conversion to succeed may turn to 
zero as the precision e turns to zero. If this is the case, we 
can not stochastically convert ps into px- The stochas- 
tic approximate convertibility seems to be the weakest 
one as far as one-copy-to-one-copy transformations are 
concerned. 

For last several years the stochastic convertibility of 



multipartite pure states has been extensively studied. In 
the case of pure states LOCC transformation reduces to 
local filtering operations The question of convert- 

ibility by local filtering operations (LFO) can be investi- 
gated using the theory of normal forms and entanglement 
monotones which are invariant under LFO, see Ref. 

Conditions for stochastic convertibility of mixed states 
are less understood. It was conjectured in Ref. [f| that a 
two-qubit mixed state ps can be stochastically converted 
into a state px iff the Bell-diagonal normal form of pT 
equals to convex sum of some separable state and the 
Bell-diagonal normal form of ps- 

In this paper we consider only mixed states of a very 
special form, namely bound entangled states associated 
with unextendible product bases in a system of three 
qubits. The notion of unextendible product basis (UPB) 
was originally introduced in Ref. @ . In our setting UPB 
is a family of vectors \Si), . . . , \S n ) € (C 2 )® 3 such that 

• Each vector \Sj) has a product form \Sj) = \Aj) ® 
\Bj) <g> \Cj) for some one-qubit states \Aj), \Bj), 
and \Cj). 

• (Si\Sj) = Sij for all i and j, 

• The orthogonal complement to the space spanned 
by the vectors . . . , \S n ) does not contain prod- 
uct vectors. 

The vectors \Sj) are refered to as members of the UPB. 
If S = {\Sj)} is a UPB, the linear space spaned by its 
members will be denoted Tis- A mixed state ps associ- 
ated with a UPB S is defined as a properly normalized 
projector onto an orthogonal complement of TLs- In our 
setting it is 
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It can not be a separable state, since its range does not 
contain product vectors. However, as was pointed out in 
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Ref. Q, in the case of three qubits ps is separable with 
respect to any bipartite cut, for instance ABC = AUBC 
(we will use letters A, B, and C to label the qubits). In 
particular it is not possible to distill some pure entangle- 
ment between any two qubits starting from many copies 
of ps and applying only LOCC transformations. For that 
reason ps is refered to as a bound entangled state. 

In section[n]we build a complete classification of three- 
qubit UPBs. It appears that each UPB has exactly four 
members. By local unitaries and permutations of the 
members we can bring any UPB into the following form: 

|5i) = |0) ® |0> <8> |0), 
15a) = 

|5 3 ) = 1^)0 11)0 1(7^), 

154) - 1^)0 1^)0 11), (2) 

where 

\A)(A\ = {1/2) [I + cos (6 A )a z + sin (9 a )(t x ], 
\B)(B\ = (l/2)[I + cos(8 B )<T z +sm(8 B )a x }, 
\C)(C\ = (1/2) [1 + cos (8 c )a z + sin (9 C )<J X ] ■ (3) 

(Here and throught the paper we use a designation \A ) 
for a state orthogonal to the state \A).) The family 
given by Eq. (|2I3|) specifies a UPB for an arbitrary triple 
9 a, 9b ,&c 7^ mod 7r. However some triples must be 
identified since the corresponding UPBs can be matched 
by local unitaries and permutation of the UPB's mem- 
bers (in this case the associated states are related by local 
unitaries). We show that the "fundamental" region of pa- 
rameters corresponds to 9a, 9b, 9c € (0,n). Restricting 
ourselves to this region we count each UPB exactly one 
time. 

Our main conclusion is that stochastic approximate 
LOCC conversion (in either direction) of states ps and 
Pt associated with UPBs S and T is impossible, unless ps 
and pt are related by local unitaries. Following Ref. 
we obtain the necessary conditions for LOCC convert- 
ibility finding the necessary conditions for convertibility 
by separable superoperators. So the statement which we 
have actually proved is following: 

Theorem 1. Suppose S and T are three-qubit UPBs 
which are not related by local unitary operators and per- 
mutation of the members. Let ps and pt be the bound 
entangled states associated with S and T. There exists 
a finite precision e(S, T) > such that for any separable 
completely positive superoperator E satisfying E(ps) =/= 
one has 

The paper is organized in the following way. In sec- 
tion |n|wc build a classification of three-qubit UPBs and 
reveal some useful facts about them. The most important 
fact is Lemma|3]which says that members of UPB are the 
only product vectors in the spanning space of a UPB. 



In section IIIII we address stochastic (exact) convertibil- 
ity of mixed states associated with UPBs. In section HVI 
we examine a simplified version of the problem, namely 
approximate convertibility by local filtering operations. 
Section |V| contains a proof of Theorem ^ In Conclu- 
sion we summarize the results obtained in the paper and 
discuss possible application of our method to UPBs in a 
system C 3 O C 3 . 

II. UNEXTENDIBLE PRODUCT BASES FOR 
THREE QUBITS 

In this section we put forward a complete classification 
of UPBs for a system of three qubits and prove some 
useful facts about them. The qubits will be refered to as 
A, B, C. Let \Sj) € C 2 O C 2 C 2 be members of UPB: 

\S j ) = \A J )®\B J )®\C j ), j = l,...,n. (4) 

By definition (Si\Sj) = Sij for all i and j. Let Hs be a 
linear n-dimensional space spanned by the basis vectors 
Sj and ps be a density operator proportional to a projec- 
tor onto orthogonal complement of Hs, see Eq. Q. As 
an example consider a UPB " Shifts" suggested in Ref. [|| . 
Its members are 

= |O)0|Q)0|O), 
\S 2 ) = |i)o|-)o|+), 
\S 3 ) = |+)o|l)o|-), 
\S 4 ) = |-)®|+)®|i). 

(Here we use a standard notations |±) = 2~ 1 / 2 (|0)±|1)).) 
If some product vector \S) = \A) \B) ® |C) ex- 
tends this basis, then one of vectors \A), \B), \C) must 
be orthogonal to at least two vectors from the family 
{|0), |1), |+), |— )}, which is impossible. 

An important characteristic of an UPB is its orthog- 
onality graph, see Ref. |8|- For our purposes it will be 
convenient to introduce more general definition of an or- 
thogonality graph. 

Definition 1. Let M = . . . ,\^ N ) £ H} be a fam- 

ily of vectors in a finite- dimensional Hilbert space Ti. An 
orthogonality graph G — (V, E) of the family M has a 
set of vertices V = {1, . . . , N} with one vertex assigned 
to each state |\E'j). A pair of vertices (i,j) is an edge, 
(ij)€E, iff(Hf i \HF j ) = 0. 

For example, an orthogonality graph of any family of 
vectors in C 2 is a collection of disjoint edges and isolated 
vertices. This definition will help us to prove the follow- 
ing statement. 

Lemma 1. Any three-qubit UPB has four members. 

Proof. First of all we note that the number of members 
n can not be greater than 5. Indeed, if n > 6, we have 
Rk(As) < 2. From Ref. @ we know that a PPT three- 
qubit state with a rank two or one must be separable. 



3 



3 3 




FIG. 1: Possible orthogonality graphs Ga for the family 



But ps is not separable by definition. Let us prove that 
an existence of a UPB with five members also leads to 
a contradiction. Let S be a such UPB with the mem- 
bers given by Eq. Consider orthogonality graphs 
G A = (V,E A ), G B = (V,E B ), and G c = (V,E C ) for the 
families {\Bj)}, and {|Cj)} respectively. Here 
V= {I,..., 5} and 

(i,j)eE A iff (A i \A j ) = 0, 
(i,j)€E B iff (B l \B 3 )=0, 
(i,j)€E c iff (Ci\Cj)=0, 

Orthogonality of the basis vectors (Si\Sj) — 5ij implies 
that any pair of vertices belongs to at least one of 
the sets E A ,E B ,Ec- Thus 

\E A \ + \E B \ + \E C \> n{n 2 1] =10 

and at least one of the sets E A , E B , Ec contains 4 or 
more edges. Assume that \E A \ > 4 and focus on the 
graph G A . Observe that G A can not have a vertex with 
a valence 3 or greater. Indeed, if, say, (1,2) S E Al 
(1,3) e E A , and (1,4) e E A then a state |S 6 ) = 
\Ai) (£> \Bi) (gi \C§) extends a basis which is impossi- 
ble. Besides, G A can not contain cycles with odd num- 
ber of edges (this constraint comes from two-dimensional 
geometry). Summarizing, G A must match the following 
restrictions: 

• there are at least 4 edges 

• a valence of any vertex may be 0,1,2 only 

• there are no odd cycles 

Up to permutations of the vertices there are only two 
graphs satisfying all the restrictions. They are shown 
on FIG. 1. For the graph on the left the basis can be 
extended by a state \S 6 ) = \A 2 ) ® \B£) <g> \Ci). For the 
graph on the right the basis can be extended by a state 
\S 6 ) = \A£)(g>\B 5 )®\C 5 ). Thus an assumption that UPB 
contains 5 members leads to a contradiction. □ 

For our purposes the order of UPB members will not 
be important. Besides we would like to identify UPBs 
which can be matched by local unitary operators. Let us 
introduce the following equivalence relation: 



Definition 2. UPBs S and S' are equivalent, S ~ S' , 
iff \Sj) = U A ®U B ® U c \S' aU} ), j = 1,2,3,4, for some 
unitary operators U A , U B , Uc and some permutation a S 

We will see later (corollary to Lemma |3J) that UPBs S 
and S' are eqivalent iff the corresponding states ps and 
ps> are related by local unitaries. Classification of UPB 
with 4 members is given by 

Lemma 2. Each equivalence class of UPBs has a repre- 
sentative of the form 

\Si) = |0> ® |0> <8 |0>, 

\S 2 ) = 

\S 3 ) = |A> ® |1> ® IC- 1 -}, 

\S A ) = \A X ) ® \B X ) ® |1). (5) 

For arbitrary \A),\B),\C) $ {|0),|1>} the set of states 
given by Eq. J3J) is a UPB. 

Proof. Suppose S = {\Sj)} is a UPB with four members 
given by Eq. J?J. The unextendability implies that nei- 
ther of the sets {|A,)}, {| .&,■)}, {|Cj-}} can contain two 
equal states. Indeed, if, say, \Ai) — \A%), then a state 
\S 5 ) = \Aj^} (8> \B^) (g) \C$r) extends the basis. One re- 
mains to take into account mutual orthogonality of the 
basis members. Applying local unitary operators we can 
make \S X ) = |0) ® |0) ® |0). It follows that each of \S 2 ), 
\Ss), \S4) contains at least one factor |1). Since neither 
of the sets {|A,)}, {\Bj)}, {\Cj)} can contain two copy 
of |1), by appropriate permutation of the members IS^), 
I1S3), |54) we can always make 

|Si> = |0)®|0)®|0), 

|S 2 ) - |1)®|B 2 >® \C 2 ), 
\S 3 ) = \A 3 )®\l)®\C 3 ), 
\Si) = |A 4 )®|S 4 )® |1), (6) 

Orthogonality (S2IS3) = translates into {C 2 \C 3 ) = 0. 
Analogously, (^l^) = implies that (B2I-B4) = 
and (S , 3 |S , 4) = gives us (A3IA4) = 0. Introducing 
\A) = \A 3 ), \B) = \B 2 ), and \C) = \C 2 ) we arrive to 
a representation declared in Eq. J^J. 

Consider an orthogonal basis {|Sj)} as in Eq. (JSJl and 
suppose that \A),\B),\C) <£ {|0),|1)}. Assume that a 
state I a) <g> |6) <S> |c) extends the basis. It is orthogonal to 
each of \Sj), j = 1, 2, 3, 4 either on the qubit A, or on the 
qubit B, or on the qubit C. But each of the states \a), \b), 
|c) may provide orthogonality to at most one member of 
the family {IS^)}. Thus the basis is unextendible. □ 

Let us count, how many real parameters we need to 
parameterize equivalence classes of UPBs. Applying local 
phase shifts to the qubits A, B, C, we can always make 

\A)(A\ = (l/2)(/ + cos(^K+sin(^K), 
\B)(B\ = (l/2)(I + cx>s(0 B y z + sin (0bK), 
\C)(C\ = (l/2)(/ + oos(0 o y + sin (0c)o")- (7) 
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Different triples (9 a, 9b, 9c) ma Y still represent the same 
equivalence class. For example, applying o~ z to the qubit 
A, we should identify 9 a and —9a (analogously, Ob = 
—9b and 9c = —9c)- Thus all equivalence classes of 
UPBs are contained in the region 9a, 9b, 9c £ (0,7r). 
Although it is not necessary for the following discussion, 
in Appendix we prove that all triples (9a, 9b, 9c) inside 
this region represent different equivalence classes. 

By definition, an orthogonal complement of a span- 
ning space Hs of any UPB S does not contain product 
vectors. The set of product vectors in the space Hs it- 
self however is not empty. It always contains at least n 
vectors — the members of UPB. Generally, the space Hs 
contains more than n product vectors (some examples are 
given in Conclusion section). Amazingly, in the system 
of three qubits the members of UPB are the only prod- 
uct vectors in Hs, which significantly simplify analysis 
of convertibility issues. 

Lemma 3. Let S be an arbitrary UPB for three qubits. 
The basis members \Sj), j = 1, 2, 3, 4 are the only product 
vectors in the spanning space Hs- 

Proof. Let S — {\Sj)} be a UPB with four members given 
by Eq. 10. We will prove more strong statement, namely, 
that the basis members \Sj), j = 1,2,3,4 are the only 
product vectors in Hs with respect to any bipartite cut 
of three qubits. Consider for example a cut ABC = 
A U BC. Suppose that for some \^a) £ C 2 and some 
I ^ bc) £ C 2 ® C 2 there exists a decomposition 

4 4 

\* A )®\* B c)=J2 a ^> En 2 = l ( 8 ) 

3=1 3=1 

Tracing out the qubits B and C we get: 

|*a)(*a| = P12 +P34, (9) 

where 

_ ( M 2 ai a 2 {B\0){C\0) \ 

912 - ^aia 2 (0|£J)(0|C> |a 2 | 2 J' 

and 

_/ |a 3 | 2 a 3 a 4 {B X \l){l\C x ) \ 

m ~ \a 3 a 4 (l\B x )(C^\l) \a 4 \ 2 J' 

Here p\ 2 is represented by its matrix elements in {|0), |1)} 
basis, while P34 is represented by its matrix elements in 
•fJA},!^!- 1 -)} basis. Since p\ 2 and ,034 are non-negative 
Hermitian operators, the equality in Eq. J§J) is possible 
only if P12 and ^34 are both proportional to the projector 
|\E , j 4)(^'a|- In particular it implies that detpi 2 = and 
det/934 = 0. According to Lemma|21the bases {\B), IB- 1 )} 
and {\C), \C )} do not coincide with the basis {|0), |1)}. 
It means that the determinant 

det Pl2 = |a 1 a 2 | 2 (l-|( J B|0)| 2 |(C|0)| 2 ) 



equals zero only if a\a 2 = 0. Analogously, det P34 = 
only if a^o^ = 0. Since the basis {\A), {A- 1 }} does not 
coincide with {|0), |1)}, the equality in Eq. J5J is possible 
only if at most one of the coefficients a±, a 2 , 03, 0:4 is 
nonzero. Thus the state \^bc) coincides with one 

of the basis members. □ 

Corollary 1. UPBs S and T are equivalent iff ps and 
Pt are related by local unitaries. 

Indeed, suppose that px = UpsU\ for some unitary 
product operator U — Ua®Ub®Uc- Then UHs = Ht 
and thus vectors U\Sj), j = 1,2,3,4 are some product 
vectors in Ht- Lemma says that these vectors are \Tj), 
up to some permutation. Thus S ~ T. The reverse 
statement is trivial. 

Corollary 2. If UPBs S and T are not equivalent, they 
can not be converted to each other by local filtering oper- 
ations. 

To prove this corollary we will need two facts pointed 
out in Ref. Q|, namely: 1) An orbit of any mixed state 
under LFO (or a closure of this orbit) contains a normal 
form, i.e. a state which has all one-particle marginals 
proportional to the identity. 2) Inside each orbit there 
is a unique (up to local unitaries) normal form. From 
the representation Eq. (JjJJ it is clear that a state ps as- 
sociated with an arbitrary UPB S automatically comes 
in a normal form. Now the statement we need follows 
from the previous corollary and the two facts mentioned 
above. To ivestigate convertibility under general LOCC 
and approximate convertibility this results however is not 
strong enough. We will need one more lemma. 

Lemma 4. Let X = Xa®Xb®Xc be a non- degenerated 
(full rank) factorized 3-qubit operator such that 

X ■ Hs — Ht 

for some UPBs S and T . Then S ~ T and X is propor- 
tional to a unitary operator: X — rU , r£l, UU^ = I. 

Proof. Obviously, multiplying X on factorized unitary 
operators on the left and on the right is equivalent to a 
choice of representatives in equivalence classes of S and 
T. So we can assume that 

• Xa, Xb, and Xc are Hermitian operators 

• Members of the UPB S are given by Eq. JSJ . 
Members of the UPB T will be denoted as 

|Ti) = |Ar) ® {Bt) ® |Ci), 

\T 2 ) = \A 2 ) ® \B 2 ) ® \C 2 ), 

\T 3 ) - ^3) ® \B 3 ) <g> |C 3 ), 

|T 4 ) = |v4 4 )®|S 4 )® |C 4 ). (10) 

If \ip) G Hs is a product vector, then its image X\tp) e 
Ht is also a product vector. Combining Lemma and 
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the fact that X is non-degenerate operator, we conclude 
that 

X\S j )=x j \T a(j) ), j = 1, 2, 3, 4, (11) 

for some coefficients Xj ^ and some permutation a G 
54. Let us concentrate on the qubit A. From Eq. (jll|l we 
infer 

X A \$) = 0x14,(1)}, 

X A \\) = 0a|4r(2)>, 
X A \A) = 3 |4r(3)>, 
X A \A^) = 0414.(4)), (12) 

where a,j are some non-zero coefficients. For each j — 
1,2,3,4 choose some unital vector |4j") orthogonal to 
\Aj). Taking into account the fact that X A is Hermitian 
and non-degenerated, we can rewrite Eq. 112(1 as 

X A \A^ 1} ) = bill), 

X A \Aj {2) ) = 6 2 |0), 

X A \A^ (3) ) = 6 3 4 X ), 

X A |4^ (4) ) - 6 4 4), (13) 

for some non-zero coefficients bj. Applying the classif- 
icaltion Lemma |2 to the UPB T we conclude that the 
family {|4j~)} and the family {|4')} coincide up to per- 
mutation of the elements and some phase adjustment for 
each element. Moreover, from this Lemma we can learn 
that an orthogonality graph for the family {|4)l mu st 
be one of three graphs shown on FIG. 2. 

am 0(2) am ao am aa) 
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a<4) a<3) CT(4) C(3) a<4) ap) 

FIG. 2: Possible orthogonality graphs for the family {|4)}- 

(We choose an ordering of vertices which depends upon 
a, since this ordering appears in Eq. (|12I13|I .') Whatever 
graph is chosen, the equations Eq. (|12I13|1 guarantee that 
an operator (X A ) 4 is diagonal both in the basis {|0), |1)} 
and the basis {4), 4^)}- By definition, these bases do 
not coincide. Thus (X A ) 4 is proportional to the iden- 
tity, (X A ) 4 — x A I, for some coefficient x A . Since X A 
is Hermitian, this is possible only if X A is proportional 
to a unitary operator. Applying the same arguments to 
the qubits B and C, we conclude that X is proportional 
to a unitary operator. Finally, Lemma |3 implies that 

□ 

III. EXACT CONVERTIBILITY OF BOUND 
ENTANGLED STATES 

Consider a pair of 3-qubit mixed states ps and px as- 
sociated with a pair of non- equivalent UPBs S and T re- 
spectively. Non-equivalence assumption implies that we 



can not convert ps into px exactly by LFO, see the corol- 
lary to Lemma |3 Now we are in position to prove more 
strong statement, namely that ps can not be converted 
into pt exactly by a separable superoperators. To sim- 
plify notations, let us denote flf a set of 3-qubit product 
operators, i.e. 

n f = {X = X A ® X B ® x c , x A , X B ,X C e L(C 2 )}. 

Lemma 5. Let {Xi £ ^/};=i,...,d be a family of product 
operators such that Yl^-i X}Xi < /. Suppose that for 
some UPBs S , T and real number p > one has the 
equality 

D 

Y^XiPsX} = PPT . (14) 
i=i 

Then the UPBs S and T are equivalent. 

Remarks: 1) The factor p can be regarded as a probability 
for conversion to succeed, since 

D 

p = Tr(Y,xjXiPs) G [0,1]. 

2) Without loss of generality we can choose D as a 
real dimension of a space of all superoperators on three 
qubits, which is D = 2 13 . Indeed, superoperators E(p) = 
XipX\ discussed in the lemma constitute a compact 
convex subset in the linear space of all superoperators on 
three qubits. 

Proof. Denote 

pi = X xPs X\. 

Since pi is a positive semidefinite operator, the equality 
Eq. ((HI) tells us that 

Bg(pO C Rg(pr), l = l,...,D, (15) 

or, equivalcntly, 

Xi ■ Rg(p,s) C Rg( PT ), (16) 

where Rg(p) is a range of the operator p. Suppose that 
for some I we have pi ^ and the operator Xi is de- 
generated. Let us show that this assumption leads to 
a contradiction. By definition, Xi has a product form, 
Xi = X A <g) Xb ® Xc, so at least one of the factors X A , 
Xb-, Xq has a rank one. Consider the case Rk(X y i) = 1, 
i.e. X A — \'4 ! ){<j>\ for some vectors 1^), |0) € C 2 . Then 

Pi = \lp){lp\® PBC, 

Pbc = X B ® X c {<t>\ Ps \<P) X B ® Xl. (17) 

Since ps is a positive partial traspose (PPT) state, the 
same does pi. Then from the equality Eq. i|17|) we infer 
that pbc is a two-qubit PPT state. But according to 
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Peres criteria [ToL Hll| it implies that p B c is separable. 
Let us choose some product vector \ip B ) \ifj c ) from its 
range. Then a product vector \ips)(S> \tpc) belongs to 
the range of pi and thus to the range of pr, see Eq. (|15p. 
Since /?t has no product vectors in its range, we conclude 
that for each I either X\ is non-degenerated or pi =0. 

Let us focus on some I with pi ^ 0. According to 
classification Lemma the spaces Rg(ps) and Kg(pT) 
have the same dimension. Since Xi is non-degenerated, 
Eq. (|TT)|l actually means that Xi ■ Rg(ps) = Rg(py), or 
equivalently, that 

X/ • H T = Ws- (18) 

According to Lemma0]this is possible only if S ~ T. □ 

Since stochastic LOCC transformations are described 
by separable superoperators, we have proved the follow- 
ing theorem: 

Theorem 2. Let S and T be unextendible product bases 
for three qubits. Consider bound entangled states ps and 
Pt associated with S and T. If some LOCC transforma- 
tion maps ps into pr with a non-zero probability, then 
the states ps and pr coincide up to local unitary trans- 
formation. 



Ms- If p = pr, it would imply that stochastic con- 
vertibility with an arbitrary small precision e is possi- 
ble (although the success probability may turn to zero as 
e — * 0). Amazingly, this is not the case. 

Lemma 6. Let S be an arbitrary UPB. Consider a set 
of achievable states Ms as in Eq. \2U\) and an arbitrary 
convergent operator sequence {Y n € Ms} n >o- Denote 

Y = lim Y n . 

n — >oo 

Then either Y S Ms or Y is separable. 

Proof. Suppose that Y ^ Ms- Obviously, in the defini- 
tion of Ms one suffices to consider normalized operators 
X. Let us agree that 

X = Xa <£> Xb <8> Xc , - \\X B \\ = \\X C \\ = 1. 

Let us choose an arbitrary sequence of normalized op- 
erators {X n S ^/} n >o, Ps[X n ] > 0, such that Y n = 
X n psXl/ps[X n ]. Without lose of generality we can as- 
sume that the sequence {X n } is also convergent (since 
X n are taken from a bounded manifold, we can always 
extract a convergent subsequence). Denote 



X = lim X n = X, 



X c 



(21) 



IV. APPROXIMATE CONVERTIBILITY: 
SIMPLIFIED SCENARIO 

Consider a three-qubit UPB S = {|Sj)} with the mem- 
bers given by Eq. J3J). Let T be another UPB which is not 
equivalent to S. In this section we will study a simpli- 
fied version of the stochastic approximate convertibility 
problem. Namely, we will consider approximate convert- 
ibility by local filtering operations, i.e. transformations 
like 



Ps 



Xp s Xl 
Ps[X] 



p s [X] = Tr(XpsX^). (19) 



where le ilf is a product operator such that ps[A] > 0. 
Our goal is to prove that there exists a finite precision e > 
such that all states p which can appear in Eq. 1191) lie 
outside e-neighborhood of pt (we use fidelity to quantify 
the distance). Denote a set of achievable states p as Ms'- 



M, 



p = ^wr ■ Xen f> ps[x}>o\. (20) 

PS [A J 



The results of the previous section imply that pr ^ 
Ms- Note however that Ms is not a compact set since 
Ps[A] = for some product operators X (for instance, 
Ps[\Sj) (Sj]} = for any member of the UPB). Thus there 
might exist a sequence of product operators {A„}„> , 
Ps[A„] > 0, such that the corresponding sequence p{X n ) 
converges to some operator p which does not belong to 



Note that ps[A] = 0, since otherwise Y S Ms- It means 
that XpsX^ = 0, or, equivalently, psX^ = which im- 
plies Im(A"t) C Ker(ps). Therefore 

lm(x\) ® Im(Ajj) ® Im^) C H s . 

But according to Lemma |3] the only product vectors in 
Hs are the members of S. Thus X = \a,b,c)(Sj\ for 
some j £ [1, 4] and some one-qubit normalized states |a), 
| b), |c) Without loss of generality, we can assume that 
j = 1, i.e. 



X A = |o>(0|, X B = \b)(0\, X o = \c)(0\, 



(22) 



(see Eq. ©). The elements of the sequence {A„}„> can 
always be written as 

X n = Xa.tl ® X B ,n ® Xc,n, 

X Atn = |a n )(0| + |a n )(l|, 
X B , n = |6„)(Q| + |/3 n )(l|, 
X C ,n = |c„)(0| + | 7 „)(l|, (23) 

where |a„), |6 n ), |c„) and \a n ), \f3 n ), |7n) are some one- 
qubit states. The requirements Eq. I|21I22() translate into 



lim \a„) = \a), lim \b n ) = \b), 



lim | c n 

n — >oo 



and 



lim (a n \a n ) = lim (/3„|/3„) = lim (jnlln) = 0. 
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Denote 



Y n = 



Ps[X n ] ' 



(24) 



To compute a limit of Y n for u->oowe can keep only the 
leading terms in the expression for X n (see a comment 
below), namely 

X w A (0) + A (1) 
a£>> = |a„,6„,c„)(0,0,0|, 

}(l,0,0| + |a„,/3„,c„)(0,l,0| 
+ K,6 B ,7 n )<0,0,l|. (25) 

Obviously, X® does not yield any contribution to Y n , so 
that 

X n psXl^X^p s X^. 

Now we CEin substitute 6^, c n in the expression for 

An by a, 6, c respectively, since we would like to keep 
only leading order terms. Using the list of the UPB's 
members Eq. (01 one can easily check that matrix 
elements (1, 0, 0|p<?|0, 1, 0), (1, 0, 0|ps|0, 0, 1) (and their 
cyclic permutations) vanish. Therefore we arrive to 

» (1,0, 0|ps]l, 0, 0) • \a n ,b,c)(a n ,b,c\ 
+<0,l,0| jOg |0,l,0)-|a,/3 n) c)(a,/3 n ,c| 
+ (0,0,l|p s |0,0,l) • \a,b, ln )(a,b, ln \. (26) 

Note that this is a separable state and that the matrix 
elements 

(1,0,0| PS |1,0,0), (0,1,01^10,1,0), (0,0,1^0,0,1), 

are strictly positive, since Lemma [3] tells us that the 
states |1,0,0), |0, 1,0), and |0,0, 1) do not belong to Us- 
It justifies that the leading terms in the expansion for 
X n psX\ do not vanish and all terms we have disregarded 
are indeed small compared with the terms we keep. In 
particular, 

Tr(xUpsX^) » 



(l,0,0|p s |l,0,0)(a n |a„) 
+<0 ) l,0|p< ? |0,l,0)(/? n |/3 n ) 
+ (0,0,l|p 5 |0,0,l)(7„| 7 „). 



We conclude that 

Y = lim Y n = lim 



Tr ( xPpsXP 1 



(27) 



(28) 



where approximations Eq. (|26I27|I should be substituted 
into nominator and denominator. But from Eq. i|26|) we 
infer that Y is separable. □ 



Denote M a closure of the set M. It is a compact set. 
A fidelity F(p, px) is a continuous function of a state p. 
According to Lemma El F(p,px) < 1 for all p £ Ms- 
Lemma implies that F(p, px) < 1 even for all p £ Ms- 
Thus there exists a finite precision e > such that 

F(p, Pt) < 1 - e, for any p £ M s . 

Of course the precision e may depend upon S and T, in 
particular e — ► as S turns to T. 



V. APPROXIMATE CONVERTIBILITY IN A 
GENERAL CASE 

We now are prepared to analyse approximate convert- 
ibility by arbitrary completely positive separable super- 
operators. Such operators correspond to probabilistic 
mixtures of LFO, so we can exploit the results of sec- 
tion llVI A proof of Theorem^is contained in the follow- 
ing lemma. 

Lemma 7. Let S and T be non- equivalent UPBs. Sup- 
pose {Xi £ f2/}i=i,...,£) is a family of product operators 
such that p s [Xi] = Tr(X lPs x}) > for all I. De- 
note ps = ^2iLiPs[Xi]- There exist a finite precision 
e = e(S,T) > depending upon S and T only, such that 



F 



<l-e(S,T). (29) 



Proof. We start from introducing an auxiliary proxim- 
ity measure between px and the final state. Let \Tj), 
j = 1,2,3,4 be the members of the UPB T. Consider a 
functional 

4 

h{p) = ^>2{Tj\p\Tj). (30) 

It is clear that 

f T (p)=0 iff Rg(p) C Rg( PT ) . (31) 

Thus the equality fx{p) — is necessary for p = px- 
Define p as 



P 



Let us first prove that there exist a finite precision 5 > 
depending upon S and T only, such that 



Mp) > s. 

Using linearity of fx we can write 



(32) 



Mp) 



Y, 



Ps[Xl, 

1=1 



ps 



MY), 



XipsX[ 
Ps [A,] 



G M s , 



(33) 



8 



see Eq. (12011 . Using the same arguments as in the 
proof of Lemma [S] one can easily show that an inclu- 
sion Rg(Yi) C Rg(p T ) is impossible for any Yi S Ms- 
From Eq. J2U we infer that fr(Yi) > for any F; € M s . 
Moreover, fr(Y) > for any Y belonging to a closure 
of the set Ms- Indeed, Lemma HO tells us that such Y ei- 
ther belongs to Ms or is separable. In the latter case the 
inclusion Rg(F) C Rg(p T ) is impossible, since Rg(pr) 
is free of product vectors. Summarizing, there exists a 
finite precision 5, such that 



fr(Y) > 5, for any Y e M s 



(34) 



Returning to Eq. I|33|l we immediately get the estimate 
Eq. 123. 

To conclude the proof we need to turn inequality 
Eq. (|32l) into an upper bound on fidelity F(px,p)- It 
can be done as follows. We start from the standard defi- 
nition of fidelity: 



f (Pt,p) = Tr (VptpVpt) 



(35) 



Introduce orthogonal projectors Pt and P^ onto the 
spanning space TIt and its orthogonal complement re- 



spectively. By definition, Pt 
(1/4) and thus 



P4 



F{pt,p) 



1 



TrJP^pP 



I. Then pt = 
I (36) 



Note that fr(p) = Tr(pPr)- From an identity 

Tr (P^pPf-) + Tr (PtpPt) = 1 
and from the estimate Eq. (|32|) we get 
Tr (P^pP^) <1-S. 



(37) 



Taking into account that the operator P^rpP^r has a rank 
at most four and applying Cauchy-Schwarz inequality we 
arrive to 



Tr ^P^rpP^ < y/4Tr (P^pP^) < 2V1 - 6. 
Substituting it to Eq. ij56*|) gives us 



(38) 



non-equivalent UPBs S and T the stochastic approxi- 
mate conversion of the associated states ps and pt is 
impossible. 

It would be interesting to apply our method to UPBs 
in some other low-dimension systems. The system of two 
qutrits C 3 ® C 3 is of particular interest, since for this 
system a complete UPBs classification has been already 
found in Ref. 0. It is known that all two-qutrit UPBs 
are characterized by the same orthogonality graphs and 
consist of five members. Unfortunately a set of product 
vectors in a spaning space of a UPB is generally larger 
than a set of the UPB's members. We have checked it 
for two particular two-qutrit UPBs, called in Ref. @ as 
" Tiles" and " Pyramid" . Using the designations of this 
reference, the extra product vectors in the spaning spaces 
of "Pyramid" and "Tiles" respectively can be written as 

|0)®|0) and i(2|0)-|l)+2|2))®(2|0)-|l) + 2|2)). 

Exploiting the symmetry of these particular UPBs one 
can show that there are exactly six product vectors in 
their spaning spaces. However the presence of extra prod- 
uct vectors as well as three-dimensional geometry makes 
the convertibility analysis very complicated. 
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APPENDIX 

In this section we prove that all UPBs given by 
Eq. with parameters 9a,Qb,8c £ (0, tt) are not 

equivalent to each other in the sense of Definition 
Consider UPBs S = {\Sj)} and S' = {\S'j)} with the 
members 



F(pt,p)<VT^5<1--- 



We have proved the lemma. 



VI. 



CONCLUSION 



(39) 
□ 



We have studied local convertibility of three-qubit 
mixed states associated with unextendible product bases. 
A complete classification of three-qubit UPBs is sug- 
gested. This family of UPBs is shown to have some nice 
mathematical properties which allow to investigate con- 
vertibility question completely. We proved that for any 



and 



|5i> 
\S 2 ) 
\Ss) 
|5 4 > 



\S[) = 
\S') = 



\s' 3 ) 
is',} 



|0)®|0)«|0), 
|1) ® \B) ® \C), 
\A) ® |1> <g> {C 1 -), 
\A X ) <g> IB 1 -) ® |1) 



|0> <8> |0> ® |0>, 
|1)®|B')®|C>. 

ll)®!^), 
\A' ± )9\B ,± )9\1). 



(A.l) 



(A.2) 
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Since the overall phase of the vectors is not important, 
the parameterization Eq. @ is equivalent to 

\A) = cos(0 A /2)|O)+sin(0 A /2)|l), 

\B) = cos(0 B /2)|O)+sm(0 B /2)|l), 

\C) = cos(0 c /2)|O)+sm(0 c /2)|l), (A.3) 

and 

\A>) = cos(^/2)|0)+sin(^/2)|l>, 

\B') = cos(^/2)|0)+sin(^/2)|l), 

\C) = cos(^/2)|0)+sm(^/2)|l). (A.4) 

Suppose that 

U A ®Ub®Uc\S 3 ) = \S a{j y) 

for some one-qubit unitary operators Ua, Ub, Uc and 
permutation cr g S 4 . If <r(l) = 1 then orthogonality im- 



plies a(j) = j for all j. It means that |(0|A)| = |(0| A')| , 
\{0\B)\ = \(0\B% and |(0|C)| = |(0|C")|. From Eq. fOjl 
and Eq. I|A.4|) we have | cos (o A /2)\ = | cos (0' A /2)\. Since 
8a,0'a £ (Oj 71 ") we conclude that 9' A = 9 a- Analogously 
one gets 0' B — Ob and 6' c = 9c- 

If a(2) = 1 then orthogonality implies er(2) = 1, 
<t(3) = 4, and ct(4) = 3. It means that \{0\A)\ = 
\{l\A' X )l \(0\B)\ = |(0|S')l> ^d |(0|C)| = |(Q|C")|. Thus 
(6 a ,6b,0c) and (6' A ,6' B ,6' C ) coincide. 

If fi(3) = 1 then orthogonality implies er(l) = 3, 
er(2) = 4, and cr(4) = 2. It means that |(0|A)| = |(0|A')|, 
|(0|B)| = 1(11^)1, and |(0|C^)| = |(0|C' X )|. Thus 
(9a,0b,O c ) and {9' A ,9' B ,9' C ) coincide. 

Finally, if er(4) = 1 then orthogonality implies = 
4, ct(2) = 3, and cr(3) = 2. It means that KOlvT 1 )! = 
|<0K X >|, |(0|^)| = |(0|B' X )|, and \(0\C)\ = |{1|C" X )|. 
Thus {9 A , 9 B , 9c) and (9' A , 9' B , 9' c ) coincide. 
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